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Abstract 



We prove the existence of solution in a class H 2 (f2) x H 1 ^) to steady compressible Oseen 
system with slip boundary conditions in a two dimensional, convex domain with the boundary of 

class H 5 / 2 . The method is to regularize a weak solution obtained via the Galerkin method. The 

0^ ■ 

problem of regularization is reduced to a problem of solvability of a certain transport equation 

by application of the Helmholtz decomposition. The method works under additional assumption 

!>■ ' on the geometry of the boundary. 
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1 Introduction 

In this paper we consider a system of Stokes-type equations describing steady flow of a barotropic, 
compressible fluid in a two dimensional, convex domain with H 5 ^ 2 - boundary, supplied with in- 
homogeneous slip boundary conditions with nonnegative friction coefficient. The system can be 
considered as a linearization of a Navier-Stokes system for compressible fluid around a constant 
flow (v = (1, 0), p = 1), thus we will call it compressible Oseen system. The slip boundary con- 
ditions involving friction enable to describe the interactions between the fluid and the boundary 
of the domain. It also turns out that they allow to extract some information on the vorticity of 
the velocity, that can be used to show that the velocity has higher regularity. Such approach has 
been applied in [5] and [7] to incompressible flows. In this paper we follow these ideas, modify- 
ing them in a way that they can be applied to the compressible system. A significant feature of 
this system is its elliptic-hyperbolic character: the momentum equation is elliptic in the velocity, 
while the continuity equation is hyperbolic in the density. Therefore we can prescribe the values 
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of the density only on the part of the boundary where the flow enters the domain and a singularity 
appears in the points where the inflow and outflow parts of the boundary meets. 

We show existence of a solution (u, w) E H 2 (Vl) x H 1 ^). A method we apply is to regular- 
ize a weak solution obtained via the Galerkin method. Analysing the vorticity of the velocity we 
can show that the density is in fact solution to a certain transport equation, obtained via elimina- 
tion of the velocity from the continuity equation. The problem of regularization is thus reduced 
to a problem of solvability of a transport equation. The values of the density are prescribed on the 
part of the boundary where the flow enters the domain, and the density can be found as a solution 
to the transport equation via method of characteristics, thus singularities appear in points where 
the inflow and outflow parts of the boundary coincides. We show that the solvability of this 
transport equation is relied with the geometry of the boundary near the singularity points, thus 
we can define classes of domains where our method of regularization can or cannot be applied. 

Since similar difficulties resulting from the mixed character of the problem appear in the 
analysis of steady compressible Navier-Stokes system, it is likely that the results of this paper 
will turn out useful in future analysis of the nonlinear problem. Now let us formulate the problem 
more precisely. 

The steady compressible Oseen system reads: 



(1.1) 



where fi is a bounded, convex domain in R 2 with a boundary T of class H 5//2 . u : il — > R 2 is 
the velocity of the fluid and w : fi — > R is the density, n denotes outward unit normal to Y. We 
assume that F E L 2 (Vt),G E H l (Vt) and B E i7 1 / 2 (r) are given functions, v and \x are viscosity 
constants satisfying v + 2/x > and / > is a friction coefficient (note that if / — > oo then 
the conditions (ll.llh d reduce to a homogeneous Dirichlet condition). The system (11.61) can be 
considered as a linearization of a steady compressible Navier-Stokes system around a constant 
flow (v = (1, 0), w = 1). More precisely, the perturbed flow satisfies inhomogeneous boundary 
conditions n ■ u\t = d and w\r in = Wi n , but if we assume that d and Wi n are regular enough 
we can reduce the problem to homogeneous boundary conditions (11.11) 4 r. Thus we distiguish 
the inflow and outflow parts of the boundary Y as the parts where the perturbed flow enters and 
leaves the domain: 

Y in = {x : m{x) < 0}, Y out = {x : m{x) > 0}. 

Let us also denote r* = {x : n\{x) = 0}. We assume that r* consist of two points: 
x* = (xi*, £2*) andx* = (x*,^) (see Fig. H). Due to the convexity of f2 we can define functions 

xi(x 2 ) and x~i{x 2 ) for x 2 E (x 2 *, x* 2 ) in the following way: 



d Xl u — /j,Au — {y + fi)Vdi 


vu + 7Vu> - 


= F in 


Q. 


divu + d Xl w = G 




in 


n, 


n ■ 2fiD(u) ■ r + /( u ■ t) = 


= B 


on 


an, 


n ■ u = 




on 


dQ, 


w = 




on 


r 

- 1 - mi 



(xi(x 2 ),X 2 ) E Tin, (xi(x 2 ),X2) E Y 



out 



Around x* and (x*) X2 is given as a if 5/,2 -function of x\. We will denote these functions by 
x l 2 {x\) and x% (xi) respectively (Fig. [2]) For convenience we will denote 
C(DATA) := C(fi, v, Q, F, G, B) The main result of this paper is 

Theorem 1. Assume that F G L 2 {Vl), G G H 1 ^), B G i/ 1 / 2 (r) and / is large enough. Assume 
further that the boundary near the singularity points satisfies the following condition 

d I < g < 3 : hm - : : = hm - : = +oo. (1.2) 

xi~^x Vt \\Xi — Xu\ q — X 2 *\ xi->s* \\Xi — X\\ q — X2I 

Then the system ( li.(5l) has a unique solution (u, w) G H 2 (Q) x H l (Q) and 

IMIfl»(n) + lk||ifi(a) < C(DATA). (1.3) 

The geometric condition (11.21) may look strange since it is formulated in a general form, but 
it has a clear meaning. Namely, the boundary near the singularity points can not be too flat, more 
precisely, our method works if the boundary is less flat than a graph of a function |xi| 9 around 
zero for some q < 3. We also show (lemma[T3~1(b)) that the method does not work if the boundary 
behaves like |a;i| 3 or is more flat. The limit case if the boundary is more flat that the graph of \x± \ q 
for all q < 3, but less flat than |a;i| 3 . An example of such a function is |a;i| 3 | In |xi||. In lemma 
[14] we show that our method doesn't work in such case. The proof of theorem [Q is divided into 
several steps. In section |2] we show existence of a weak solution in a class (H 1 ^)) 2 x L 2 (f2) 
using the Galerkin method (theorem [3]). To obtain a weak solution it is enough to assume that 
G G L 2 (Q) , and no further constraint on the geometry of T is required. The constraint (11.21) arises 
when we want to show that the weak solution belongs to class H 2 {VL) x H l (Vt), and we also need 
G G H l {Q). The issue of regularity of the weak solution is treated in section First we prove 
that the vorticity of the velocity belongs to H l {Vt) (lemmafT]). Such approach has been applied to 
incompressible Navier-Stokes equations in [5 J and [7J. In the incompressible case we can next 
solve a div-rot system to show higher regularity of the velocity, but in the compressible case we 
have to extract some information on the density. The idea is to use the Helmholtz decomposition 
in H l {Q), that means express the velocity u as a sum of a divergence-free vector function and 
a gradient. The standard theory of elliptic equations enables us to show that the divergence-free 
part belongs to H 2 (tt), and in order to show higher regularity of the gradient part it is enough to 
show that div u G i/ 1 (fi). In lemma ITOl we show that div u + w G H 1 (fi), thus we have to show 
that w G H 1 (fi). The method is to show that the density is a solution to the transport equation 

jw + w Xl =H G H\tt). (1.4) 

Thus the problem of regularization of the weak solution is reduced a problem of solvability of 
the transport equation (11.41) . The boundary condition (11.61k prescribes the values of the density 



on the inflow part of the boundary and (11.41) can be solved via method of characteristics, thus 
a singularity appears in the points x* and x*, which we will call the singularity points. It turns 
out that we can solve the equation (11.41 ) provided that the singularity is not too strong, what is 
reflected in the constraint (1 1 .2ft . We will finish the introductory part removing inhomogeneity on 
the boundary. Let us construct a function u G H 2 (tt) satisfying 

n-2/iD(w) -t + f(u-r)\ r = B and n ■ u\ r = 0, (1.5) 



such that ||«o||i/ 2 (n) < C(^)||-6||ifi/2(r)- Then a pair [u, w), where u = u — u , satisfies 

d Xl u — fiAu — (u + ^Vdivu + 'y'Vw = F in Q, 

div u + d xi w — G in Q, 

n ■ 2/xD(u) • r + /( u ■ r) = on T, 

ra • tt = on T, 

w = on r»n, 



where 



F = F + /^Amo + {y + /j,)Vdivu - d Xl u e L 2 {9) 
G = G-divu eH^Q). 



(1.6) 



(1.7) 



Obviously we have 

II^IU^) < C(\\F\\ L2m + \\B\\ L 2 (r) ) and \\G\\ mm < C(\\G\\ mm + ||B|| L 2(r)), 
thus from now on we can work with the system (11.61) denoting u := u, F := F , and G = G. 



(X 1 *,x 2 *) 




Out 



(X^*,X2*) 



Figure 1 : The domain 



2 Weak solution 



In order to define a weak solution to the system (11.61) consider a space 

V = {v e C°°(n) : v ■ n\ r = 0, n ■ 2fiD(v) ■ r + f(v • r)| r = 0} 



X 2 -x 2 (x.|) 




Figure 2: The domain, functions xi, Xi, x%, x\ 



and V — Vq Hl<n) equipped with the norm | \v\\v 



vWh 1 ^)- Consider also a space 



W = { V e L 2 (n) : Vxi G L 2 (n) and r 1 \ Tm = 0} 

with the norm | \w\\ w = \\w\\ L2{n) + \\w Xl \\ L2{n) . 

Now we want to introduce a weak formulation of (11.6b . First, observe that for u, v regular 
enough we have 

/ (—fiAu — {y + /i) V divu) ■ vdx = / 2/jD(m) : V v + v div u div v dx — 
Jn Jn 

/ n ■ [2/iD(w)] ■ v da — I n- [u(divu)Id] ■ vda, (2.1) 

where A : B = £^ ;=1 a y fey for A = {oy}, B = {a hl } e R nxn . 

Thus taking u eVq'vs\ (11.6b i and multiplying it by a function v G Vo we get 

/ {f • 9 Xl u + 2/iD(w) : Vf + v divu div v — ^w divv} dx + f(u ■ r) (v • r) da = 
Jn Jv 

= F-vdx. (2.2) 
Jn 

Multiplying (11 .6b ? by a regular function 77 G W we get 

?7[dw « + w xi ] dx — / Gr/dx. (2.3) 



The above considerations leads to a natural definition of a weak solution to the system (11.6b . 



Definition 1. By a weak solution to the system di.(5D we mean a couple (u,w) G VxW satisfying 
dZ2D - tiZMfor each (v, rf) E V x W. 

We want to show existence of a weak solution using the Galerkin method. In order to show 
existence of solutions to approximate problems in section l2TT1 we apply well known result (lemma 
[[])• This result automatically gives uniform boundedness of the sequence of approximate solu- 
tions, what enables us to show convergence of approximate solutions to the weak solution in 
section [2T2l 

2.1 Approximate solutions 

In order to construct a Galerkin approximation of a weak solution let us introduce an orthonormal 
basis ofV: {4>i}°^i = {(4>j, (ff)}^ and finite dimensional spaces: V N = {J2i=i a i4>i '■ a i £ Tt} C V. 
We will search for a sequence of approximations to the velocity in the form 

N 
.N 



E C ^- (2-4) 

Let us denote x\ := Xi(x 2 ). Taking u = u N , v = <p k and w = w N where 



u 

i=l 



w N (xi,x 2 )= / (G — divu N )(s,x 2 ) ds 



in (|2.2I) we get 



V c? I d xi <j>i ■ <P k dx + 2/i V tf / D(&) : V<p k 
i Jn i Jn 

+v\^ cf / div <fii ■ div <fik dx — 7 2_, { {G ~ /, c i div<pi)(s,X2)ds} div<j)kdx 
i Jn i Jn Jan i 

+f Y, cf [(fa ■ t) (0 fc • r) da = I F- <\> k dx.(2.5) 
i Jt Jn 

For k = 1 . . . N we obtain a system of N equations on coefficients {cf}f =l . If a function u N of 
a form (12.41 ) satisfies the equations (12.51 ) for k = 1 . . . N, it means that a pair (u N , w N ) satisfies 
(|2.2|) - (|2.3I) for each (v, 77) e V^ x W. We will call such a pair (m^, w*) an approximate solution 
to d2j) - (1231) . 

The system (|2.5I) . fc = 1 ... AT is rather complicated thus in order to solve it we will use the 
following well known result (see for example (9l|): 

Lemma 1. Let X be a finitely dimensional Hilbert space and let P : X — > X be a continuous 
operator satisfying 

3M>0: (P(£),£)>0 for ||£|| = M (2.6) 

Then 3£* : ||£*|| < M and P(f*) = 



In order to apply lemma [Q we will need some auxiliary results in spaces V and W. 
Lemma 2. (Poincare inequality in V) 

VveV: \\u\\ LH n) < C{n)\\Vu\\ LH n). (2.7) 

Assume that (12.121) doesn't hold. Then 3{v k }^ =1 G V such that 

Proof. WV Vk\\ L 2(p\ < j \\vk\\L 2 (n)- Without loss of generality we can assume ||ffc||L 2 (n) = lV/e, 

thus 

||Vv fc ||r»(O)->0. (2.8) 

Clearly {t> fc } is a bounded sequence in if x (f2) and thus thanks to boundedness of VL the compact 
embedding theorem implies that it contains a subsequence {v k } that is a Cauchy sequence in 
L 2 (f2). But (12.81 ) implies that Vv k . is also a Cauchy sequence in L 2 (H). Thus {t> k .} is a Cauchy 

sequence in if 1 (f2), hence v^ — ► u* for some w* G iJ 1 (fi). Obviously ||u*||L 2 (n) — 1 an d 
||Vv*|| = 0, thus v* is constant almost everywhere. But also (t>* ■ n)\r = 0, and since O is a 
bounded set with regular boundary, the unit normal takes all the values from the unit sphere on 
T. Therefore 

a.e. \ 

v = const _^ * a - e - n 
<■ => v = 0, 
(v*-n)|r = 0j 

what contradicts 1 1 v*\\ L 2 ^ = 1 D 

Now we will use the Poincare inequality to show that in V a following modification of the 
Korn inequality holds: 

Lemma 3. Assume that f is large enough. Then for u G V : 

f 2/iD 2 (w) + I f(u-r) 2 da> C\\u\\ 2 m . (2.9) 

Jq Jt 

Proof. The proof is based on a proof of a different version of the Korn inequality in [5 ]. We have 

2 „ k 



2 [b 2 ( U ) = J2 [«) 2 +<<] = iiv«ni» (n) + [ y £ui.ui i dx= 

= 1 1 Vu| ||2 (n) + J2 U ^ u k dx- J2 u * < ni da - / XX ui ra i da - (2 - 10) 
^ij=l • /r i,i=l ' Jr i,j=l 

The second term of the r.h.s is equal to J n div 2 udx > and the third term vanishes since 
(u - n) |r = 0, thus from (12.101) we get 

k 

2a [ B 2 (u) > fi\\Vu\\ 2 LHn) -ti [ f" v^u* ni t da, (2.11) 

J si J? i „-_i 



but we have | J r ^ . =1 w* m- 7 n^. dcr| < C(fi) | |w| \l 2 (t) an d thus using the Poincare inequality 
(1277b we get 

/ 2 D 2 (u) + /(« ■ r) 2 > C(fi, /i) | H |iri(n) + [/ - C(fi, fi)] \\u\\ LHr) 
Jn 

and the last term will be positive provided that / is large enough. □ 

The last inequality we need is the Poincare inequality in W. 
Lemma 4. (Poincare inequality in W) 

Vr]eW: \\T]\\ L2 (n) < diam(Q)\\r] Xl \\ L2{n) . (2.12) 

Proof. The proof is straightforward using density of smooth functions in W and the Jensen 
inequality. □ 

The following theorem gives a solution to the system (12.51) 

Theorem 2. For F,Gg L 2 (Q) and B e L 2 (T) there exists a solution {cf}^ to the system 
d2.5D . k — 1 . . . N. The function u N = J2i c f 4>i satisfies 

\\u N \\ m{n) <C(DATA). (2.13) 

Proof. In order to apply Lemma Q] we have to define an appropriate operator 

p n . v n ^ v n_ For convenience let us den ne B N :V N xV N ^ R: 

B N (Z N , V N )= f v N d Xl i N + ^ f D(£ N ):Vv N + v [ dtv£ N divv N 
Jn Jn Jn 

-1 I {f \G-divi N ){s,x 2 )ds}divv N dx + f I \i N ■ r)(v N ■ r) da - f F ■ v N dx. 
Jn Jxx_ Jr Jn 

Now (12.51) can be rewritten as B[u N , 4>k) = and thus it is natural to define 

P N (Z N ) = J2 BN ^ N ^^ for ^6V W . (2.14) 

i 

We have to verify the assumptions of Lemma[TJ Obviously P N : V N — > V N and it is a continu- 
ous operator. For £ = £\ af^i we have 

N N 

(p N (e),n = czb^^m^, j>f^) = 

fc=l i=l 

AT TV JV 

= Ei 5 ^ &) E a ? (**> **) > = E 5iV ^ **)°? = 5A ^ *")■ ( 2 - 15 ^ 



fe=l i=l fe=l 



Using the definition of B N we can rewrite (12.151) as 

(P N (Z N ),S N ) = 2/i / D 2 (£ N ) dx + v [ div 2 i N dx+ f £ N d xl £ N dx + f f(£ N ■ r) 2 da 
Jn Jn Jn Jv 



h 



"7 /{ / \G-dwi N ){s,x 2 )ds}dw£ N dx- I F-f dx. 
Jn Jx± Jn 



h 



Using the Korn inequality (12.91 ) we get I1 + I2 > C(fi, u, fi) 1 1^1 \ 2 Hl for / large enough. Now let 
us denote 

r] N (xi,X2)= / (G — div ^ N )(s,x 2 ) ds. 

J X\ 

Then 77^ = G — div £ N and we have 

J 3 = - 7 I ' r, N divi N dx= [ v N vZ~ [ G V N dx> 
Jn Jn Jn 

> f G V N dx > -C\\G\\ L 2 (n) \\ V N \\ L 2 {n) > -C\\G\\ L2{n) (\\G\\ L2{n) + H^Hi^ft)). 
Jn 

Combining these bounds we get 

(P N (Z N ),n > C(ji,n)\\t N \\ 2 m ia) - (\\F\\l^) + ||G||r"(0))||eil^(n) - \\G\\l Hn) . 

Thus there exists C = C(ji, fi, F, G) such that (P N (£ N ), £ N ) > for 1 |f 1 1 = C, and apply- 
ing lemma [Q we conclude that 3£* : P N (£*) = and ||£*|| < C. But since {(pi} is a basis 
of V, the definition of P N (12.141) yields 

P N (C) = «=* B N (t, 4> k ) = 0,k = l...N, 
thus £* is a solution to (1231) D 

2.2 Existence of weak solution 

Now we show that the sequence (£ N , r) N ) constructed in previous section converges to the weak 
solution of our problem. 

Theorem 3. Assume that F,Ge L 2 (Q) and f is large enough. Then there exists a weak solution 
(u,w) to A1.6\) satisfying the estimate 

\\u\\ v + \\w\\w<G{DATA). (2.16) 



Proof. The estimate (12.131) together with (12.121) gives 



\u 



N\\ , IL.JVll , ||„,JV 



hhci) + IKIIl^) + ||<||^(a) < C(DATA). (2.17) 



Since the sequence {w x } is bounded in L 2 (f2), there exists a subsequence u;^* and a function 
( G L 2 such that w^ fe — ^ £. Now let us denote for simplicity w N := w k . It is bounded in L 2 , 

t-2 

thus there exists a subsequence w Nj — ^ iu for some function w E L 2 . Now we need to show that 
C = w xi , but this is quite obvious. We have 

\/v e L 2 : - w Nh v Xl = / w^ fc u -> / C v and / /» « Xl - 
Jn Jn Jn Jn 

thus f n ( v = — j n w v Xl Vf G L 2 (fl) . 

It is a bit more complicated to show the existence of u. The estimate (12.171) gives boundedness 
in L 2 (f2) of the sequences {div u N } , {d Xl u N } , {D(u N )} and boundedness in L 2 (T) of {(u N ■ 
r)}. Thus up to a subsequence 

divu N ££, d xl u N ^a, D(u N )^(3 and u " . r ^ 5 
for some £,a,(3 G L 2 (fi) and some 5 G £ 2 (T). On the other hand, since the sequence {u N } 

t-2 

is bounded in H 1 , the compactness theorem yields u N — > u up to a subsequence for some 
u G L 2 (Vt). We want to show that in fact u G H l and that (u, w) satisfies (12.21) - (|2.3I) . But we 
have V</>GC^: 

- / ud Xl <f>< / u n d Xl (f)= / (f)d Xl u n ^ a<f> 

Jn Jn Jn Jn 

thus a = d Xl u. Similarily we can verify that 

£ = div u 

(3 = T>{u) (2.18) 

5 = u ■ t\t- 

Thus it G H\Q), and the pair («, w) satisfies (Q - (Q ViV G N V(v, 77) G V* x W. The 
density of V^ in 1/ implies that it also satisfies (12.21) - (12.31) V(t> , i))gVx W 7 . Thus indeed (it, w) 
is a weak solution. The estimate (12.161) is obtained in a standard way taking v = u and rj — w 
in (|2.2I) - (12.31) and then applying the Korn inequality (12.91) and the Poincare inequality in W 
(12121) . □ 

3 Regularity 

In this section we will show that the weak solution belongs to a class H 2 (Vl) x H 1 ^). The idea 
of the proof has been outlined in the introduction. We start with showing that if (u, w) is a weak 
solution then rot u G H l {Vl). 

10 



Since on this level we have only weak solutions, we have to work with the weak formulation 
(12.21) - (|2.3I) . Consider a special class of test functions: 

V 1 = {veV: V^: (j>eH 2 (n), u • n| r = 0, 0|r = 0} 

where V -1 = (d X2 , —d Xl ). Note that on T we have ■£ — v • n — 0. 



2 -u 1 

xi res 

Since div v = for «gVl, thus for t> G VI (12.21) takes the form 



Let us denote a = rot u = u — u x 



[ ad Xl (pdx + 2fi I D(u) : Vvdx = [ F-V ± <pdx- f f{u-r)-^-da (3.1) 
Jn Jn Jn Jv on 



Lemma 5. For u G V, v G V 1 we have 



/ 2uD(u) : Vf dx = — fi / aA(f)dx+ / 2 (/ix — /)(«• r)—— da (3.2) 

Jn Jn Jv on 

where a = rot u and x denotes the curvature of T. 

To prove lemma[5]we will use following auxiliary result, proved in 0: 
Lemma 6. For u G Vo we have 

rotu\ r = (2 X --)(u-T), (3.3) 

ft 

where x is the curvature ofT. 

Proof of lemma\5\ Due to density of Vo in V it is enough to proove (13.21) for 
u t G Vq, v G V\. For such functions we have (we omit the subscript e): 



2/iD(-u) :Wvdx = — j 2fidivD(u) ■ v dx + / n ■ 2/fD(u) ■ v da 

Jn Jv 

Since we have 2divD(u) = A u + Vdiv u, and using the definition of Vo we can write 

2fjbD(u) : Vvdx = - / fi{Au + V divu) ■ V ± <pdx- f (u ■ r) -£- da. (3.4) 

n Jn Jv on 

Integration by parts yields 

/ Vdivu ■ \7 ± (j)dx = / divu—-da = (3.5) 

Jn Jv or 

and 

/ Au ■ V^(j) = 4> Arot udx+ (pAu-rda 
Jn Jn Jr 

V «, ' 

=0 

V0-Vrotw+ / (f) —-rotu = / rotuAcj) — I rot u-— da (3.6) 

o Jt on J n J T dn 
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Substituting (1331) and d3T6b into (|3~4l) we get 

/ 2/iD(w) : Vfdx = — /i \ rotu A(f)dx + // rotu— -da — / f (u ■ t) da, (3.7) 
in in ir on i r 

and application of (13.31) to the boundary term yields (13.21) □ 
With lemmal5l (13.1l) takes the form 

ad Xl <f)dx- aA<pdx= F-V ± <pdx- [ (2lix - f)(u ■ r) -£- da. (3.8) 
in Jq in ir un 

Since u G H 1 ^), we can construct d G H 1 ^) such that 

d|r= (2/*x -/)(«• r ). (39) 

IMIIff^n) < CIMIh 1 ^)- 

Now consider a decomposition a = b + d where b\r = 0. From (|3.8I) we see that the function b 
satisfies 

bd Xl (f)dx+ / Vb-V(f)dx = - / dd Xl <pdx+ I F-V X (j)dx- I Vd-V<f>dx. (3.10) 
in in in in 

Inverting the above reasoning we can prove 

Lemma 7. Assume that (u, w) is a weak solution to ( li.(5D . Then rotu G H x (f2) and 

||roH| ffl(n) <C(/JATA). 

Proo/ Consider a problem: find 6 G -f^(fi) satisfying (13.101) V0 G -f^(fi). Obviously this 
problem has a solution 6 G Hq(Q) satisfying 

\\b\\ H i < C(DATA,d) < C(DATA,\\u\\ H i) <C(DATA). 

In particular b satisfies (13.101) V0 G ifg H i? 2 . Thus if we define a* = b + d, where d is given 
by ([3T9i then ||a*||fli < C(.DA7M) and a* satisfies (ES) V0 G Fq 1 n H 2 . But this means that 

a* = rot u □ 

We will use this fact together with a well known result, the Helmholtz decomposition in 
H l {tt) (©,©): 

Lemma 8. (Helmholtz Decomposition) For u G if x (f2), fnere exxste fie H 2 (Q) such that 
n ■ A |r = and 

M = V^ + V X A (3.11) 

Now our goal is to show that if (w, u>) is a solution to (12.21) - (12.31) then ■0, A G H 3 (tt), thus 
u G /7 2 (ft). 

Lemma 9. Assume that (u, w) is a weak solution to ( |i.(5|) and (■0, A) is the Helmholtz decompo- 
sition ofu. Then A G H z (Vt) and \\A\\ H : i{n) < C(DATA). 
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Proof. On the boundary we have n ■ A 1 - = r ■ \7A = fjA, thus the condition n ■ A 1 - \qq = yields 
A\dQ = const. Moreover, rotu = rot(X7ip + A- 1 ) = rot A 1 - = AA. We see that A is a solution 
to the following boundary value problem: 

AA = a 
A\r = const, 

where a = rot u E H 1 (fi). Since the boundary is of class H 5 ^ 2 , the standard theory of elliptic 
equations yields A E H 3 (£l) and | |-A]]ir 3 (r2) < C(^) Hulls' 1 to) '-' 

Now we want to show that also ip E H 3 (tt). We have divu = Aip and on the boundary we 
have = u ■ n = (Vip + V ± A) • n = Vip ■ n. Thus ip satisfies 

Aip = div u 

^lr = 

and in order to prove that ip E H 3 (tt) it is enough to show that div u E H 1 (fl). The next step is 
to prove the following 

Lemma 10. Assume that (u,w) is a weak solution to the system rti.<5D . Then 

- (2fi + v) div u + -f w =: H E H\Q) (3.12) 

and 

|| - (2/i + v)div u + -f w\\ H i {n) < C(DATA). (3.13) 



Proof. For u E Vq we can integrate by parts in (|2.2I) and using (12.11) we obtain 

/ {v ■ d Xl u — \p Au + [y + [i)Vdivu] ■ v } dx — / / ywdivvdx = / F-vdx. (3.14) 
Jn Jn Jn 

Substituting the Helmholtz decomposition to (13.141) we get 

/ — [y + 2//)V(A'0) • v dx — \ r )wdivvdx = 
Jn Jn 

[ (F - d Xl (Vip + V ± A) + fiAV^A) -v dx. (3.15) 

Jn^ v ' 



F 



From lemma|9]we see that F E L 2 (f2) and ||F|| L 2(q) < C(DATA). Integrating formally by 
parts the second term of the r.h.s. of (13.151) we get 

\-(2{i + v)V{A4))+-iVw\-vdx= F-vdx. (3.16) 

Jn Jn 

At the beginning we assumed that u E Vq in order to write (13.141) . but we can understand the 
identity 2D(w) = Am + Vdiv u in a weak sense, and thus we have 

F = V[-(2// + u)Aip + yw] = V[-(2/i + v)divu + w] 

and so lemma ITOl is proved □ 
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Combining (13.121) and (11.61) 9 we get 



W + Wxi = H +G=:HeH\Q). (3.17) 



v + 2fj, x 2fi + v 

We see that the density is a solution to a transport equation. Our goal is now to use this fact to 
show that w G iJ 1 (fi). Since we already know that w G L 2 {VL), the problem reduces to showing 
that w X2 G L 2 {Vt). A natural way to extract some information on w X2 from the equation (13.171) is 
to differentiate it with respect to x 2 . For simplicity we will write 7 := -^— and H := H. We 
have 7W + w Xl = e~ yxi d Xl {e ,xx w), thus differentiating (13.171) with respect to x 2 we get 

d X2 [e-^d Xl (e^w)} = d X2 H G L 2 (Q). (3.18) 

We want to use the above identity to define w X2 in an appropriate way. In order to do this assume 
first that w X2 G L 2 {Vt) is well defined. Then (13 . 1 81) can be rewritten as e~ 1Xl d^ \e 1Xl w„ 2 ] = d X2 H, 
thus 

d Xl [e^w X2 ] =e^d X2 H=:a (3.19) 

If we assume also that w X2 is well defined on T in , then we can write: 

e^ Xl w X2 (x u x 2 ) = e' y ^ X2) w X2 (x 1 (x 2 ))+ f a{s,x 2 )ds. (3.20) 

J Xl(x2) 

This identity will enable us to define w X2 on tt provided that it is well defined on T in . The 
boundary condition w\r m = implies that the tangent derivative of w is well defined on T in : 
-§pw\r in = 0. Provided that the first order derivatives of w are well defined on T in , this identity 
can be rewritten as 

T lr w Xl + t 2 w X2 = 0, (3.21) 

but due to (13.171) we have w xi \r in = (w + w Xl )\r in = H\ Tin G H 1/2 (T in ), and thus (13.211) can be 
rewritten as w X2 \r in = — ^ H \r in . Note that on r^ we have ^ = xi_(x 2 ) (Fig. |3]), thus we can 
rewrite (T3.20I) as 

w X2 (x 1 ,x 2 )=e~ JXl [-e lx ^ {x2) Xi(x 2 )H(x 1 (x 2 ),x 2 )+ / a(s,x 2 )ds\. 

Jxxixi) 

Since a G L 2 (Vt), we see that 

pxi, rxl(x 2 ) 

w X2 (xi,x 2 ) G L 2 {Vi) <^=^> / / [xi'{x 2 ) H(xi(x 2 ),x 2 )] 2 dxi dx 2 < 00. 

Jx 2 * J x\(xi) 

Since the function under integration doesn't depend on X\, we can rewrite the r.h.s. of the above 

as 

r x 2 

/ [xi{x 2 ) - xi(x 2 )][xi(x 2 )] 2 H 2 (xi(x 2 ),x 2 ) dx 2 . 

J X2* 

For simplicity let us denote H(x 2 ) := H(xi(x 2 ),x 2 ), The above considerations leads to the 
following conclusion: 
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Xi(x 2 ) 



Figure 3: The function £1(0:2) and the tangent vector 



Lemma 11. Let (u, w) be a weak solution to A1.6\) . Assume that 

r x 2 

/ (3(x 2 )\xi(x 2 )\H 2 (x 2 )dx 2 < 00, 

•J X2* 

where 
and define 



(5(x 2 ) = [xi(x 2 ) - xi{x 2 )}\xi(x 2 )\ 



X(xi,x 2 ) 



7^1 r _ pTEiC^)™ ' 



x 1 (x 2 )H(x 1 (x 2 )) + 



.''1 



a[s 



xi(a;2) 



x 2 ) ds\ . 



(3.22) 



(3.23) 



(3.24) 



where a is defined in A3.19\) . 

Then X € L 2 {Vt), ||A|| L 2 (n) < C(DATA) and X = w X2 . 

Remark. The functions H(xi(-)) and J* 1 ,* a(s, ■) ds are defined a.e. in (x 2 *,x 2 ), thus A is 
defined a.e. in Cl, more precisely, it is defined for all x\ and almost all x 2 E (x 2 *, x^). 

Proof of lemma [771 Since a e L 2 {VL), we see that (13.221) implies A G L 2 {VL). Moreover, 
inverting the passage from (13.181) to (13.201) we conclude that 

e-^d Xl [e' rXl X} = d X2 H, 

thus indeed A = w X2 . D 

Now we are ready to formulate a regularity result that can be considered a major step in the 
proof of theorem [TJ 
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Proposition 1. Let (u,w) G V x W be a weak solution to di.(5D and assume that the boundary 
constraint t\3.22\l holds. Then (u,w) G H 2 (Q) x H l (Q) and 

IN|fl»(ft) + lklk(n) < C(£>ATA). (3.25) 

Proof. At this stage in order to complete the proof it is enough to resume the steps we have 
made. From (|3.17l) and lemma [TT] we have w G H 1 ^). Thus from (|3.12l) we conclude that 
divu G H l (VL), and so (13.121) yields ip G H 3 (Q), where u = V^ + V^A From lemma [9] we 
have A G H 3 (Q), hence we conclude that u G H 2 (Q) and the estimate (13.251) holds. □ 

As we see, the condition (13.221) is crucial for our regularization method to work, but it is 
hard to interprete it since it doesn't depend only on the geometry of the boundary, but also on 
the function H. Thus we want to formulate some conditions equivalent, or at least sufficient for 
(13.221) to be satisfied, that would depend only on the geometry of T. Such condition is stated in 
the following 

Lemma 12. Assume that for some e > 

/ 2 p 1+e (x 2 )\x[(x 2 )\dx 2 < ex), (3.26) 

then (TJ[22]) holds. 

Proof. Since the integrability in (13.221) is questionable only in the neighbourhood of x 2 * and x 2 , 
we can fix some small S > and focus on 

CX2*+S 



rX2*+o px 2 

( / + / )f3(x 2 )\x^(x 2 )\H 2 (x 2 )dx 2 . 

J X2* J x%—5 



We will consider the first integral, the second is dealt with in the same way. Observe that on T in 
we have dx 2 = \r 2 \ da, thus 

/ l3{x 2 )\x{{x 2 )\H 2 {x 2 )dx 2 = [ ' (3{x 2 )\^\H 2 {x 2 )dx 2 ~ [ p H 2 da, 

Jx2, JX2* T JT} n 

where T\ n denotes the part of T in between x 2 * and x 2 * + o. In the last passage we used the fact 
that n — 1 in the neighbourhood of the singularity points. Since H G H 1 / 2 (r in ), due to the 
Sobolev imbedding theorem we have H G L p (T in ) Vp < +oo, and thus 



3e > : / /3 i+e da < oo 
rL 



(3H Z da <+oo, (3.27) 



but on T] n we have da ~ \xi_{x 2 ) \ dx 2 and the l.h.s of (13.271) is equivalent to 

rX2*+& 



3e > : / (3 1+e {x 2 )\xi(x 2 )\dx 2 < +oo 

J X2* 



u 
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The condition (13.261) depends only on the geometry of T in the neighbourhood of the singu- 
larity points. Now we want to determine some classes of domains where the condition (13.261) 
holds and doesn't hold. We will focus on one of the singularity points, let's say x 2 * and assume 
without loss of generality that (xu,x 2 *) = (0, 0). For simplicity let us denote x l 2 (xi) =: l(xi). 

To start with, consider a class of domains where l(xi) = \xi\ q , q > 2. We have to assume 
q > 2 to assure that \x\ q E H 1 / 2 (R). Indeed, we have 

\x\ g E tf^R) *=* [ [ l( x + ^~ xr l 2 dhdx< +00, (3.28) 

Jo Jo h 

where r = q — 2. Dividing the integral over x into f Q + J h we can see that r.h.s of (13.281 ) is 
equivalent to integrability on (0, 1) of a function x 2r "" 1 , what holds for r > 0. 
We have xi(x 2 ) = —x 2 and xT(x 2 ) = x 2 , and thus 



3+2£ eg 1 A 
(3 1+e (x 2 )\x 2 ;(x 2 )\ = {xT(x 2 ) - x 1 (x 2 )] 1+e |xi'(x 2 )| 2+6 ~ x 2 " 



We see that (I3T261) holds for 



g < ^±^£ < 3 (3.29) 

H 1 + e V 

In particular, (|3.26l) doesn't hold for any e > (or even for e = 0) if q — 3, but for any q < 3 

there exists e such that (13.261) is satisfied. Although this example concerns only a particular class 

of boundaries, it suggests that we should be able to determine whether (13.261 ) holds or does not 

hold by comparing the function l(xx) with the limit case from our example, i.e. l*(xi) = \xi\ 3 . 

Let us denote 

g q =\im l p±. (3.30) 

xi -»0 \Xi\ q 

It turns out that whether (13.261) holds depends on g q in the following way: 

Lemma 13. Let g q be defined in t\3.30\l . Then we have 

(a) 3q < 3 : g q = +00 => (\3.26\) holds for some e > 0; 

(b) #3 < +00 =>- ( 13.2(51) does not hold for any e > 0. (3.31) 

Proof. Let us show (b). We have 

1 1 ,- i dxAxi) i r | ^ 2 (4 /3 ) 1 ,. l d X2 {x 1 2 /3 ) 1 

thus (we understand that ~ = 00 and — = 0): 







From (|3.32l) we get 



lim |xi'(x 2 ) I = -j — r lim \x 2 |. (3.33) 

x 2 ^0+ |<73| X2-+0+ 



1/3 1/3 

\g 3 \= lim I 2 |= lim J 2 |, (3.34) 

x 2 ^0+ Xi{X 2 ) x 2 ^0+ Xi{X 2 ) 
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thus 

1/3 



lim \xi(x 2 ) — Xi(x 2 )] = -. — r lim x 2 ■ (3.35) 

rr 2 — 0+ | (73 1 x 2 ^0+ 



Combining (13.331 ) and (13.351 ) we get 



U+^U l(™ M 2 + e 



lim [xi(x 2 ) - xj_(x 2 )] \xi(x 2 )\' 

22—0+ 



1 



-- — — — lim \xl /3 \ 1+t \x 2 2/3 \ 2+e = -. — rr— - lim \x 2 \ x 5. 

1^313+26^^0+1 * I 12 1 |^3 | 3+2e 12-O+ ' 

what implies J [xi(x 2 ) — Xi(x 2 )] 1+<E |xi / (x2)| 2+<E (ix2 = +oo since -j — Jj+ar > 0. Thus (b) is 
proved, (a) can be shown exactly in the same way by comparing l(x\) with the function \x\ \ q . D 

A remaining question is what happens in the limit case when g q = V 1 < g < 3 but 
g 3 = +oo. The following lemma gives the answer: 



Lemma 14. 



V X < <? < 3 9q ° 1 =► B16\) does not hold for any e > (3.36) 

93 — +oo J 



Proof. First of all, observe that 



x 2 \ 1/q ,. |x 2 | 1/9 



|g g | = lim - — - — — = lim 



x 2 ^ 0+ |Xi (X2) I ^2—0+ |xi(x2)| 

.V 3 

For a given function Xi(x 2 ) let us define a function h as |xi(x 2 )| = ^rr. Then we have 

1/3 
+ 00 = lim - — - — — = lim h(x 2 ) (3.37) 



22— 0+ Xi(X 2 ) X2— 0+ 



and 



We have 



x 



1/9 



3- 9 



V 1 < q < 3 : = lim 1 — 2 —- = lim x 2 39 h(x 2 ). (3.38) 

22— 0+ pi(»2)| £'2—0+ 



x 2 2 ' 3 h(x 2 ) - ti(x 2 )xl /3 ,x 2 2 ' 3 [h(x 2 ) -x 2 h'(x 2 )} 

l ^ {X2)l ~ ' ^53 = ' ^5 ' 

thus for e > 0: 






[h{x 2 )-x 2 h'{x 2 )f h 

V 

s e (x 2 ) 




18 



In order to determine whether this function is integrable when x 2 — ► 0, observe first that A e (-) is 
not integrable. Indeed, (13.381) implies x 2 h 5+3e (x 2 ) — •> V r > 0, thus for x 2 small enough 



-1-- —l— - 

2 ^ u '2 _ ~ 1_ 



and if we choose r < ~ the last function is not integrable, and thus A e (-) is not integrable. Now 
let us see what happens with B € (-). We have 

lim x 2 h\x 2 ) = lim #^-= lim -^M 
x 2 ^o+ 22-^0+ (m(x 2 )) x 2 ^o+ m(o; 2 ) 

thus /i(x 2 ) is dominating in -B(a; 2 ) when a; 2 — > 0, and in particular lim, E2 _^ + B e (x2) = +oo. We 
conclude that 

Ve>0: / A e (x 2 ) B € (x 2 ) dx 2 = +oo 
Jo 

what completes the proof. □ 

Lemma fl4l together with point (b) from lemma [T3l shows that if 
g = V 1 < g < 3 then (3 <£ L 1+€ (T in ) for any e > 0, thus we can not show (13.221) without ad- 
ditional information on the function H; the only information that we have under the assumptions 
of theoremEQis that H E H l ' 2 (T m ). 

The condition from the point (a) of lemma [T3l means that the singularity in Xi{x 2 ) in the 
neighbourhood of the singularity points cannot be too strong, more precisely, it must be weaker 
than the singularity of d X2 (x 2 ) for some q < 3. In other words, the boundary around the 
singularity points cannot be too flat, it must be "less flat" that a graph of a function \x\ \ q for some 
q < 3 (after an obvious translation). Examples of domains that allows or does not allow the 
application of our method are shown in Fig. HI The proof of our main result is almost complete. 

Proof of theoremUl If31<g<3: g q = +oo, then from lemma IT3l (a) we see that (13.261) 
is satisfied, thus lemma (fT2l) gives (13.221) . and so proposition [Q yields that the weak solution 
(u, w) G H 2 (n) x H 1 (fi). We want to show that (u, w) satisfies (11.61) almost everywhere. 

Clearly (12.31) implies that (11.61) 9 is satisfied a.e. Taking a test function v E V x H^ (VL) we see 
that also (|1.6l) i holds. The definition of spaces V and W implies that boundary conditions <\1.6\) a 
and (|1.6I) r hold, thus it is enough to show that also fll.6l h is satisfied. Since u E H 2 (Q), we can 
integrate by parts the r.h.s of (|2.2I) and obtain \/v E V: 

[F — {d Xl u — ii Au — [y + /i) V div u)~\ ■ v dx = / [n ■ 2/iD(n) ■ r + f(u ■ r)] (v ■ r) da, 
Jn Jt 



S v- 

=0 



thus indeed n ■ 2/iD(m) ■ r + f{u ■ r) a = 0. 

We have shown that for F E L 2 (Q) and G E H 1 ^) the system (11.61 ) has a solution 
(u, w) E H 2 {Q) x H 1 ^). Now let u E H 2 be and extension of the boundary data (11.51 ) and 
let (u, w) be a solution to (11.61 ) with F = F and G = G defined in (11.7b . Then (m + m , w) is a 
solution to (ll.lh and the estimate (|1.3h holds. D 
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3q<3:g q =+oo 

g 3 <+ °° 




Figure 4: Behaviour of the boundary near (x u , x 2 *) 



4 Conclusions 

We have shown existence of a solution (u, w) G H 2 (Q) x iJ 1 (fi) to the compressible Oseen sys- 
tem with slip boundary conditions (11.11) . The method we applied follows the approach of 0, Q 
and reduces the problem of regularization of the weak solution to a problem of solvability of the 
transport equation (11.41) . We can solve this equation and thus prove that the density w E H 1 (il) 
provided that the boundary constraint (11.21 ) holds. It should be underlined that this constraint 
does not result from the system (11.61) itself, but from the method of regularization that reduces 
the problem to solvability of (11.41) . Application of different methods of regularization might en- 



able us to weaken the assumption (11.21) . In particular it would be interesting if we could weaken 
it in the way that enables domains where ri\ = on a set of positive measure, where clearly 
(11.21) cannot hold. A natural continuation of this paper would be to consider the compressible 
Navier-Stokes system. A similar approach enables again to reduce the problem of regularization 
of the weak solution to solvability of a transport equation, which is however more complicated 
than (|1.4p since it contains a nonlinear term u ■ Vu>. A possible way to solve this equation is to 
apply a method of elliptic regularization. 

We also plan to extend the approach presented in this paper to L p - framework. 
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